Spin-polarized electric currents in quantum transport 
through tubular two-dimensional electron gases 
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Scattering theory is employed to derive a Landauer-type formula for the spin and the charge 
currents, through a finite region where spin-orbit interactions are effective. It is shown that the 
transmission matrix yields the spatial direction and the magnitude of the spin polarization. This 
formula is used to study the currents through a tubular two-dimensional electron gas. In this 
cylindrical geometry, which may be realized in experiment, the transverse conduction channels are 
not mixed (provided that the spin-orbit coupling is uniform). It is then found that for modest 
boundary scattering, each step in the quantized conductance is split into two, and the new steps 
have a non-zero spin conductance, with the spin polarization perpendicular to the direction of the 
current . 

PACS numbers: 72.25.-b,72.20.Dp,72.25.Mk 
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I. INTRODUCTION 

The conductance of a confined system connected to 
electronic reservoirs (held at slightly different chemical 
potentials) can be expressed by its scattering properties, 
notably the transmission, leading to the celebrated Lan- 
dauer formula^— When the motion of the electrons is 
ballistic and the transmission is perfect, the two-terminal 
conductance is given by 
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(1) 



where N± is the number of conducting channels not in- 
cluding spin. As a function of the energy of the injected 
electrons, additional channels open successively and the 
conductance follows a staircase structure, as was indeed 
observed in experiment^ When the motion is not bal- 
listic, Eq. ([T]) is modified, with the number of channels 
being replaced by the sum over the transmissions among 
channels 

Here we extend this picture to include spin-orbit inter- 
actions, and derive both charge and spin conductances in 
terms of the scattering matrix of the mcsoscopic system. 
Our approach is based on scattering theory in which the 
entire system is represented by its scattering matrix^ 
We avoid the ambiguities associated with the definition 
of spin currents in systems with spin-orbit interactions^ 
by considering scattering through a finite region where 
these interactions arc effective, and computing the cur- 
rents far away from that region, where these interactions 
are absent^ In this respect our approach differs from 
those presented in Refs. H and [h], which have deduced 
the conductance from the energy spectrum, or of Ref. UTI . 
which employed semiclassical arguments to study spin 
polarized currents in smooth barriers. 

When the electrons are restricted to move in a plane, 



an asymmetry in the confining potential leads to the 
appearance of the Rashba spin-orbit interaction^ An- 
other source for the spin-orbit coupling is the Dresselhaus 
mechanism^ in particular the cubic one, which might 
be quite substantial in GaAs nano- structures.—^ The 
Rashba spin-orbit interaction is in particular interest- 
ing, since its strength can be controlled in experiment, 
by applying gate voltagesJ^— From the strictly theo- 
retical point of view, however, there is not much dif- 
ference between the linear Dresselhaus interaction and 
the Rashba term, as they are connected by a unitary 
transformation^ The scattering formalism presented in 
Sec. [H]is independent of the type of the spin-orbit cou- 
pling. 

Spin-orbit interactions, which couple the momentum 
of the electron to its spin, have attracted much interest 
due to the possibility to manipulate the spin by electric 
fields. When electrons are injected from spin-polarized 
electrodes into a spin-orbit coupled region, they might 
lose their polarization; this decohcrence effect has been 
studied in Ref. The role of lateral interfaces formed 
between regions with different spin-orbit couplings 2 ^ in 
polarizing the spins, and its anal ogy with optics, has 
been established in Refs. [22] and [22J. However, as di- 
rect measurements of the electron spin polarization or 
of the electron spin accumulation are not easily accom- 
plished (see however Ref. |24| ) , it may be helpful to study 
the effect of that polarization on the (more easily ac- 
cessed) electronic charge conductance. At the same time, 
it is desirable to analyze the spin currents which arise 
upon the application of a small source-drain bias and ex- 
plore their possible effect on the charge transport. Pre- 
vious studies of the conductance through a region where 
spin-orbit interactions are active had relied on numeri- 
cal approximations j 2 ^— and consequently some of the 
conclusions drawn were dictated by a specific choice of 
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parameters. 

Here we study coherent transport through a mcso- 
scopic hollow cylinder; this geometry enables an exact 
calculation of the scattering matrix, and reveals its sym- 
metries. In particular we are able to analyze the ef- 
fect of these symmetries on the spin polarization of the 
transport current. We present in Sec. HI] our scatter- 
ing formalism for spin-dependent transport which follows 
closely the derivations of the Landauer formula for spin- 
independent potentials as formulated in Rcfs. [5] and H 
We identify there the specific combinations of the ele- 
ments of the scattering matrix that convey the informa- 
tion on the spin-polarization direction and its magnitude. 
By carrying out this analysis we are able to show that the 
transmission matrix of the finite system gives directly the 
direction and the amount of the spin polarization. More- 
over, under not too stringent conditions the conductance 
itself carries information on the spin polarization. In Sec. 
IHII we study the implications of this generalized Landauer 
formula for a specific system: a coherent mesoscopic hol- 
low cylinder. This system, whose experimental realiza- 
tion is found in the carbon nanotubes,— has recently at- 
tracted much interest since the remarkable observation 
of the effect of the electron motion on the direction of its 
spin.— Another experimental realization of this geome- 
try may be found in the core-shell nanowiresi^ Finally, 
we summarize our conclusions in Sec. IIVI 

II. LANDAUER FORMULA FOR 
SPIN-DEPENDENT TRANSPORT 

In scattering theory^ the electron field operator (a 
spinor) is expressed in terms of the scattering states of 
energy E, 

*(r,t) = / ^ ,H Ec M ,(£) Xm ,(r;£) . (2) 

ann 

(Here and below h = 1.) In Eq. J2]), x a nri( r "iE) i s a 
solution of the Schrodinger equation 

(E-H) Xann (r;E) = 0, (3) 

with % being the Hamiltonian of the entire system. It 
is the scattering state excited by an electron of spin po- 
larization r\ incoming from the nth channel of terminal 
a. The operator c anr] (E) \c\ nrj {E)] destroys (creates) an 
electron in such a state. The thermal average of these 
operators (denoted by a bar) is defined by the tempera- 
ture T and the chemical potential of the reservoirs. For 
reservoirs of un-polarized electrons, 



cL v (E)c alnW (E') = 2*6(E - E')5 aa ,5 nn ,5 m ,f a (E) , 

(4) 

where 

f a (E) = (e( E -^ k ° T + l)- 1 (5) 



is the Fermi distribution with the chemical potential fi a 
of reservoir a £2- Equations ([2]) and ((4]) are used to obtain 
the thermal averages of the charge and the spin currents. 
This is accomplished by writing down the currents in 
terms of the field Eq. and then performing the 
thermal average over the creation and destruction oper- 
ators, c and c' , according to Eq. 

To avoid cumbersome notations, we confine ourselves 
to the case of two terminals, located for concreteness 
along the x direction. Then, all currents flow along x 
and in the absence of a magnetic field are given by 

anr} 

X ^((Xan 77 ( r ;- E )k J IH^)Xa„ J) (r;-E))+Cc) , (6) 

where cr- is the jth Pauli matrix. The charge current, 
I x , is given by Eq. (J6j) upon identifying (Tj—a as the unit 
matrix (and multiplying by the electronic charge e). For 
a two-terminal system, the lead index is a = L or R, for 
the left and for the right reservoirs, respectively. 

The currents given by Eq. ([6]) are computed far away 
from the scattering region, where they are uniform, e.g., 
on the left side. To this end we represent the scattering 
states by the incoming and outgoing exciting waves, e.g., 

Xan V {x in L, y, z; E) = 5 0)L w^(r; E) 

+ J2 UJ °Ln'r,'(r;E)S Ln , v , tanri (E) . (7) 

n'rj' 

Here, w^° ut are the exciting waves, incoming or outgoing 
in channel n of lead a and having the spin polarization 77, 
normalized to carry a unit flux, and S is the scattering 
matrix of our system. The meaning of Eq. ((7|) is quite 
transparent: the scattering state which has been excited 
by a wave incoming in lead a consists, when considered 
in the left lead L, of all waves that are scattered into 
that lead from lead a (including those that are reflected) , 
as represented by the second term in Eq. (|7|), and the 
incoming wave in a, in case a coincides with L.— 

Next, the expansion ([7]) is inserted into Eq. ©, and 
the integration over the cross section is performed. As a 
result, I J X becomes spatially-independent, 

a nn' 

= [^(fLW-fR(E))T*'E l M Ln ,, an (E)vj, (8) 

nn' 

where the trace is carried out in spin space. Here we have 
introduced the definition 

M Ln ',an(E) = 5 Ln ,, an (£)$L,a« ( E ) - ( 9 ) 

where SLn'.Rn denotes an entry of the scattering matrix 
which is a 2x2 matrix in the spin space. The second 
equality in Eq. ((5J) is derived from the unitarity of the 
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scattering matrix. Note that I J X vanishes when there is no 
bias voltage, /i L — fi R . Although it appears as if pertains 
for general values of the bias voltage, Eq. ([8} is strictly 
valid only in the linear-response regime, in where fiL — Hr 
approaches zero. A detailed account of the subtleties 
associated with the Landauer formula can be found in 
Ref. 0. 

The general expression Eq. ([8]) for the charge and the 
spin currents is our central result in this section. In the 
small-bias limit it yields the linear charge and spin con- 
ductances. To express those, it is useful to note that the 
transmission from channel n to channel n' is determined 
by the matrix A4, Eq. @. Since this matrix is obvi- 
ously hermitian, it can be decomposed into a scalar and 
a vector^ 

M Ln ,, Rn {E) = l(r n , n (E)+V n , n (E) ■ a) , (10) 

where both T n 'n and V„/ n are real. The two real eigen- 
values of this matrix are 

*n'n( E ) = \{r n , n {E) ± \V n , n (E)\) . (11) 

Their corresponding eigenvectors are spinors fully polar- 
ized along ±V n '„. The total transmission between these 
two channels is given by 

Tn'n = Tr[M Ln , tRn ] = A+„ + \~, n , (12) 

while the spin transmission is 

Tr[M Ln ,, Rn <T] = V„, n , (13) 

with 

|V„g = A+ n -A;, n . (14) 

As is well known, (see e.g., Refs. [Ill and I22I) spin-orbit 
interactions may turn one of the modes (for each en- 
ergy) to be evanescent; when this happens, the eigen- 
value belonging to that mode, A - , is vanishingly small. 
As a result, the (perfect) conductance is reduced to 
A + ~ 1 (in units of e 2 /h), while the current becomes al- 
most fully spin-polarized, its respective conductance be- 
ing |V„„'| ~ A„ n , ~ 1. On the other hand, when none 
of the modes is evanescent, one has A„ n , ~ X~ n , ~ 1, 
leading to Q = 2e 2 /h and zero spin polarization for the 
perfect conductor. 



III. CHARGE AND SPIN TRANSPORT 
THROUGH A COHERENT TUBULAR 
TWO-DIMENSIONAL ELECTRON GAS 

In order to illustrate the general features of the scat- 
tering matrix as embodied in Eqs. (fTU)) - ([l"4")) . we study 
a specific example: transport through a two-dimensional 
stripe, of width 2d (along x), in which the electrons are 



subject to the Rashba spin-orbit interaction. The inter- 
faces at x = ±d are parallel to the y— axis; effects of 
interfacial scattering are included by a repulsive delta- 
function potential located at the interface. Elsewhere, 
the electrons move ballistically. Such a potential is char- 
acterized by a single parameter, £ (measured in momen- 
tum units). When £ is very large, the interface ap- 
proaches the tunnel-junction limit. This model becomes 
particularly transparent when a periodic boundary con- 
dition along the y axis is assumed; then the spin-orbit 
interaction does not mix the transversal modes. This 
model system, whose realization may be found in nan- 
otubes (see in particular Ref. [13) is depicted in Fig. [T] 
For cylinders, the coordinate y represents the azimuthal 
angle around the cylinder (with perimeter C), and the 
coordinate z is radial, i.e. perpendicular to the surface 
of the cylinder. 

a=0 a?0 a=0 

( ( ( 

1 1 > 

-d d x 

FIG. 1: (Color online) A schematic representation of our 
model: a hollow cylinder, containing a region, \x\ < d, where 
a spin-orbit interaction, of strength a, is active. The elec- 
trons are moving on the surface of the tube ballistically, and 
thus mode mixing is hindered. Short-range boundary scatter- 
ing, of strength separate the spin-orbit coupled region from 
the free regions. The boundaries in-between these regions are 
marked by thick solid lines. 

The scattering matrix is computed by solving for the 
scattering states of the Hamiltonian 

H = ^— + -(s(x + d)+ 6(x - dj) 

+ ^- (®(d ~ x)Q(x + d)p ■ (z x er) + he) , (15) 
2m V / 

where p is the (two-dimensional) momentum operator 
and a is the spin-orbit coupling (in momentum units). 
With the periodic boundary condition, of period C, the 
wave functions vary along y as exp[igy], and the trans- 
verse momentum is quantized, 

q = ±2nn/C , n = 0, 1, 2, . . . . (16) 

Below, we measure all momenta (including a and £) in 
units of h/ C. The x component of the wave vector is 

k = yj2mE - q 2 , \x\ > d , (17) 

for both spin components in the absence of the spin-orbit 
interaction, and 

fc„ = j(p s ±a?-q* , (18) 
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for waves in the range \x\ < d where the spin-orbit inter- 
action is effectived Here, 



Ps 



= \J 2mE + a 2 . 



(19) 



As q 2 < 2mE, the wave vector k u is always real; in con- 
trast, ki is purely imaginary for \q\ > p$ — a (we assume 
a > 0), and then one of the waves in the region \x\ < d 
is evanescent] 11 ' 22 

Following the standard procedure of matching bound- 
ary conditions at the two interfaces x = id (allowing for 
the delta-function potentials) 21 ' 33 the scattering matrix 
of each channel takes the form 



S = -1 
where 



(S + -S_)a x 
a x (S + + S_)a x 



(20) 



k + 2iC + iX± - i(T z (q - Z±) + i<r x Y ± 



(21) 



The functions X±, Y±, and Z± are all real, independent 
of whether fc; [sec Eq. (|18[) ] is real or imaginary, 



A 4 



Y, 



Ps ( s A C u(Ps ~ oQ ±g] + SujCiiPs + a ) Tg] 
Ps ( s u [<lCi T (Ps + «)] s i b°u ± (Ps - «)] ' 



Ps_ ( 



Ps 



where 



ki k, 



=1 - c,c. 



pg — a 2 + q 2 
k t k u 



S l S u 



Here we use the shorthand notations 

Si,u = sin(2fc, u d) , C Lu = cos(2fc, u d) 



(22) 



(23) 



(24) 



When the transversal channels are not mixed, then the 
total charge and spin transmissions are given by 



T = 2J T nn i and V = ^ V rm , 



(25) 



see Eqs. ©, (fT2]l. and (H31). The scattering matrix Eq. 
([20]) is unitary— and is self-dual 3 ^ upon changing q to 
—q. (Other symmetries of the scattering matrix, in par- 
ticular in the presence of the Zeeman interaction, are dis- 
cussed in Refs. [36| and H3-) The self-duailty symmetry 
of the scattering matrix implies that upon summing over 
all channels, as indicated by Eqs. (|2"5)) [i.e., over positive 
and negative values of the transverse wave vector q, see 
Eq. (fT6))] the contributions of the x and z directions of 
the polarization are cancelled, and the net spin transmis- 
sion V is along the y direction, normal to the direction 
of the current. 11 ' 28 



Figures [2HH were drawn using Eq. (|2T)1) . The curves 
in Fig. [2] are computed for entirely transparent bound- 
aries (£ = 0); indeed, the transmission in the absence 
of the spin-orbit coupling has the familiar perfect stair- 
case structure, Eq. |T|). Due to the choice of periodic 
boundary conditions in the transverse direction, the per- 
fect transmission without the spin-orbit interaction be- 
comes 2N±_ = 2(2n + 1), where n = 0, 1, 2, ... is 
the channel number. The transmission computed in the 
presence of that coupling follows roughly this pattern, 
but with two distinct new features. Firstly, it shows in- 
terference oscillations. These arise since the spin-orbit 
interaction plays the role of a potential step at \x\ < dM- 
One notes the complete vanishing of the spin conduc- 
tance in the leftmost part of the plot. There, the energy 
is too low to support a nonzero transverse momentum 
q, and consequently the motion becomes effectively one- 
dimensional. In such a case, one may handle the effect 
of the spin-orbit interaction by a gauge transformation 
which multiplies the wave function by exp[2ipsd). Since 
this phase factor is the sole effect of the spin-orbit inter- 
action, the spin-polarization vanishes complctelyi 3 ^ (Os- 
cillations as a function of the spin-orbit coupling strength 
have been discussed in Ref. [40-) When q ^ and the 
motion restores its two-dimensional character, this gauge 
transformation is no longer possible, and indeed there 
appears spin conductance. The oscillations in the charge 
conductance, resulting from the potential-step aspect of 
the spin-orbit interaction, do persist. 
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FIG. 2: (Color online) The transmission T through a stripe 
with spin orbit interaction (thick line), the y— component 
of the spin transmission V y (dotted line), and the staircase 
structure obtained in the absence of the spin-orbit coupling 
(thin line) in a perfect system, as a function of the energy. 
Here, the spin-orbit strength is a — 0.9n. 

The second prominent feature that the spin-orbit in- 
teraction induces in the conductance pattern originates 
from the possibility to have evanescent waves. At the be- 
ginning of each of the n ^ steps [see Eq. (fTE]) ]. one of 
the waves pertaining to that energy becomes evanescent, 
since using Eq. (TT5)) one has 



in 



(26) 
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where k is real. Then, C; and — iSi of Eqs. (|24[) are 
both of order 0.5 exp[2«c(|, making the determinant of 
the transmission matrix A4, Eqs. (fT0|) and (|20|) . expo- 
nentially small. This implies in turn that one of the eigen- 
values, Eqs. (|TT|) . i.e. A~ n , is very small. As a result, it 
follows from Eq. (|12[) that the (perfect) transmission is 
reduced to T nn ~ A+ rl ~ 1. Concomitantly, the polariza- 
tion magnitude, |V„„(g)|, [see Eq. HU] is also equal to 
A+„ ~ 1. However, upon adding V(q) + V(—q) the re- 
sulting vector V is along the y— axis, as shown in Fig. [5] 
This y— component decreases gradually from towards 
— 1 as the energy increases in the steps where T — An. 
At higher energies, both k u and ki are real, and the (per- 
fect) transmission becomes T nn = A+ rl + A~„ « 2, while 
the spin transmission becomes |V„„(g)| = A+ n — A~ n 0. 
For the parameters used in Figs. [5] and the contribu- 
tions of all the steps except the highest one are in this 
latter regime, and therefore each of them contributes 2 
(if q = 0) or 4 (if q ^ 0) to the transmission and to 
the spin transmission. Thus, the charge conductance ex- 
hibits 2(2n + 1) and An multiples of e 2 /h, while the spin 
transmission is roughly between and — 1. For higher 
values of a, one of the waves may remain evanescent un- 
til the next step begins, yielding only An steps in the 
transmission.— In our dimensionless units, this condition 
becomes a > (2n + l)n/n. 
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tance is not quantized, the spin-polarization is small (but 
still negative) for all values of the energy. 
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FIG. 4: (Color online) Same as Fig. [2l but with stronger 
boundary scattering, £ = 4, and a — .9n. 



IV. CONCLUDING REMARKS 

In conclusion, we have derived a Landaucr formula 
for the charge and the spin conductances. The corre- 
lations between the two conductances and their quan- 
tization have been demonstrated for transport through 
cylinders with stripes which have nonzero spin-orbit in- 
teractions. Measurements of the charge transmission can 
thus yield information on the spin polarization. It would 
be interesting to chock these predictions experimentally, 
in particular on the new set-ups of core-shell nanowires 
or made of carbon nanotubes. 
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FIG. 3: (Color online) Same as Fig. [2] but in the presence of 
potential barriers at the boundaries, £ = 1.2. Here, a — 0.9n. 

Repeating the computation for finite barriers at the 
two interfaces [see Fig. [3] shows that a modest amount 
of interface scattering is not detrimental, and the two 
main features discussed above are still detected. Even 
more interesting are the curves shown in Fig. @] Here, 
the amount of scattering at the interfaces has been in- 
creased such that the staircase structure of the charge 
conductance, both in the presence and in the absence 
of the spin-orbit coupling, is smeared; since the conduc- 



Acknowledgments 

We thank P. Silvestrov for helpful comments. AA 
and OEW acknowledge support by the binational Israel- 
United States foundation, by a grant from the Ger- 
man Federal Ministry of Eduction and Research (BMBF) 
within the framework of the German-Israeli Project Co- 
operation (DIP), and the hospitality of NTT Basic Re- 
search Laboratory at Atsugi-shi, Japan. YA acknowl- 
edges partial support from the ISF, and the hospitality 
of the Graduate School of Mathematical Science at the 
University of Tokyo. 



* Also at Tel Aviv University, Tel Aviv 69978, Israel 



1 R. Landauer, IBM J. Res. Dev. 1, 233 (1957); Philos. Mag. 



6 



21, 863 (1970). 

2 Y. Imry and R. Landauer, Rev. Mod. Phys. 71, S306 
(1999). 

3 Y. Imry, Introduction to Mesoscopic Physics, 2nd ed. (Ox- 
ford University Press, Oxford, 2002). 

4 D. A. Wharam, T. J. Thornton, R. Newbury, M. Pepper, 
H. Ahmed, J. E. F. Frost, D. G. Hasko, D. C. Peacock, D. 

A. Ritchie, and G. A. C. Jones, J. Phys. C21, L209 (1988); 

B. J. van Wees, H. van Houten, C. W. J. Beenakker, J. G. 
Williamsom, L. P. Kouwenhoven, D. van der Marel, and 

C. T. Foxon, Phys. Rev. Lett. 60, 848 (1988). 

5 M. Biittiker, Phys. Rev. B 46, 12485 (1992). 

6 Y. Levinson, Phys. Rev. B 61, 4748 (2000). 

7 J. Shi, P. Zhang, D. Xio, and Q. Niu, Phys. Rev. Lett. 96, 
076604 (2006). 

8 B. Nikolic, L. P. Zarbo, and S. Souma, Phys. Rev. B 73, 
075303 (2006). 

9 A. V. Moroz and C. H. W. Barnes, Phys. Rev. B 60, 14272 
(1999). 

10 Y. V. Pershin, J. A. Nesteroff, and V. Privman, Phys. Rev. 
B 69, 121306(R) (2004). 

11 P. G. Silvestrov and E. G. Mischenko, Phys. Rev. B 74, 
165301 (2006). 

12 E. I. Rashba, Fiz. Tverd. Tela (Leningrad) 2, 1224 (1960) 
[Sov. Phys. Solid State 2, 1109 (I960)]; Yu. A. Bychkov 
and E. I. Rashba, Pis'ma Zh. Eksp. Teor. Fiz. 39, 66 (1984) 
[Sov. Phys. JETP Lett. 39, 78 (1984)]. 

13 G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

14 J. J. Krich and B. I. Halperin, Phys. Rev. Lett. 98, 226802 
(2007). 

15 A. A. Kovalev, M. F. Borunda, T. Jungwirth, L. W. 
Molenkamp, and J. Sinova, Phys. Rev. B 76, 125307 
(2007). 

16 T. Koga, J. Nitta, T. Akazaki, and H. Takayanagi, Phys. 
Rev. Lett. 89, 046801 (2002). 

17 M. Konig, A. Tschetschekin, E. M. Hankiewicz, J. Sinova, 
V. Hock, V. Daumer, M. Schfifer, C. R. Becker, H. Buh- 
mann, and L. W. Molenkamp, Phys. Rev. Lett. 96, 076804 
(2006). 

18 T. Bergsten, T. Kobayashi, Y. Sekine, and J. Nitta, Phys. 
Rev. Lett. 97, 196803 (2006). 

19 E. I. Rashba and V. I. Sheka, in Landau Level Spectroscopy, 
G. Landwehr and E. I. Rashba, eds. (Elsevier, Amsterdam, 
1991). 

20 B. Nikolic and S. Souma, Phys. Rev. B 71, 195328 (2005); 
R. L. Dragomirova and B. K. Nikolic, Phys. Rev. B 75, 
085328 (2007). 

21 V. A. Sablikov, A. A. Sukhanov, and Y. Ya. Tkach, Phys. 
Rev. B 78, 153302 (2008). 

22 M. Khodas, A. Shekhter, and A. M. Finkel'stein, Phys. 
Rev. Lett. 92, 086602 (2004); A. Shekhter, M. Khodas, 



and A. M. Finkel'stein, Phys. Rev. B 71, 125114 (2005). 

23 V. Marigliano Ramaglia, D. Bercioux, V. Cataudella, G. 
De Filippis, and C. A. Perroni, J. Phys.: Condens. Matter 
16, 9143 (2004). 

24 I. Adagideli, G. E. W. Bauer, and B. I. Halperin, Phys. 
Rev. Lett. 97, 256601 (2006). 

25 F. Mireles and G. Kirczenow, Phys. Rev. B 64, 024426 
(2001). 

26 Shi-Liang Zhu, Z. D. Wang, and Lian Hu, J. Appl. Phys. 
91, 6545 (2002). 

27 M. Governale and U. Ziilicke, Solid State Comm. 131, 581 

(2004) . 

28 M. Eto, T. Hayashi, and Y. Kurotani, J. Phys. Soc. Jpn. 
74, 1934 (2005). 

29 The spin-orbit interaction in nanotubes is not predomi- 
nantly of the Rashba type used in Sec. Mil see T. Ando, 
J. Phys. Soc. Jpn. 69, 1757 (2000). However, the results of 
Sec. [H] are valid for these systems. 

30 F. Kuemmeth, S. Hani, D. C. Ralph, and P. L. McEuen, 
Nature 452, 448 (2008). 

31 L. J. Lauhon, M. S. Gudiksen, D. Wang, and C. M. Lieber 
Nature 420, 57 (2002); J. Noborisaka, J. Motohisa, S 
Hara, and T. Fukui, Appl. Phys. Lett. 87, 093109 (2005) 
G. Zhang, K. Tateno, T. Sogawa, and H. Nakano, Appl 
Phys. Express 1, 063003 (2008); J. W. W. van Tilburg, R 
E. Algra, W. G. G. Immink, M. Verheijen, E. P. A. M 
Bakkers, and L. P. Kouwenhoven, Semicond. Sci. Technol 
25, 024011 (2010). 

32 The extension of the formalism presented in Sec. [II] to the 
case of spin-polarized leads is rather straightforward. 

33 When writing down the boundary conditions, one should 
match the covariant momentum-which includes the spin- 
orbit contribution. 

34 Since S^ 1 + (5^) _1 = 2, the scattering matrix Eq. is 
unitary. 

35 C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

36 F. Zhai and H. Q. Xu, Phys. Rev. Lett. 94, 246601 (2005). 

37 P. Brusheim and H. Q. Xu, Phys. Rev. B 74, 233306 
(2006). 

38 This effect resembles in a way the appearance of Andreev 
reflections at the interface between a superconductor and 
a normal metal, see G. E. Blonder, M. Tinkham, and T. 
M. Klapwijk, Phys. Rev. B 25, 4515 (1982). 

39 A. A. Kiselev and K. W. Kim, Phys. Rev. B 71, 153315 

(2005) ]; note that Sanchez and Serra [Phys. Rev. B 74, 
153313 (2006)] do not find polarization for an harmonic 
confining potential along y. 

40 A. G. Mal'shukov, V. V. Shlyapin, and K. A. Chao, Phys. 
Rev. B 66, 081311 (2002); C-H Chang, A. G. Mal'shukov, 
and K. A. Chao, Phys. Lett. A 326, 436 (2004). 



